* RMatch1975y2001control.sps.
* Match patients with teeth filled in March 2000 against other
* teeth treated in year ending March 2001.
* Written by PSKL on 28/11/01.
* Amended on 17/02/03.

* Restrict March 2000 data to fillings.

get file = 'D:\PSKL\Rteeth197510003to10103.sav'.

select if ( fsched = 10003).

select if ( any( trtcode, 1401,1402,1403,1404,1421,1426)).

* save master file of filled teeth (temp1).

Compute Fillrec=1.

save outfile='c:\temp1.sav'.

get file='c:\temp1.sav'.

* save look-up table, one record for each tooth, for matching (temp2).

Compute Filled=1.

sort cases by surname initial sex dobirth tooth quadrant doacc docomp.

Aggregate outfile=* /presorted
 /break surname initial sex dobirth tooth quadrant
 /filled =max(filled)
 /dofacc dofcomp=first(doacc docomp).

save outfile= 'c:\temp2.sav'
 / keep surname initial sex dobirth tooth quadrant filled dofacc dofcomp.

* Create lookup file of patients who have had fillings in March 2000 (temp2A).
* This file contains a dummy quadrant to indicate that at least one tooth in
* the tooth position has been filled.

Aggregate outfile=*/presorted
 /break surname initial sex dobirth tooth
 /dofacc dofcomp=first(dofacc dofcomp).
compute doacc=dofacc.
compute docomp=dofcomp.

string quadrant (A2).
compute quadrant='AA'.

compute fillrec=-1.

save outfile='c:\temp2a.sav'.


* Match against quadrant records for next year, but exclude correction records.

get file = 'D:\PSKL\Rteeth197510003to10103.sav'.

sort cases by surname initial sex dobirth tooth quadrant doacc.

Match files file=*
 /table 'c:\temp2a.sav'
 /by surname initial sex dobirth tooth.

Execute.

select if (not(missing(dofacc)) and not(dofacc=doacc)).

Compute fillrec=0.

Execute.

* Remove quadrants which had fillings in March 2000.

Match files file=*
 /table='c:\temp2.sav'
 /by surname initial sex dobirth tooth quadrant.

If (missing(filled)) filled=0.
select if (filled=0).

Execute.

* Combine with March 2000 dummy quadrant records, when patients received fillings
 in at least one quadrant, to get dates of March 2000 courses of treatment.

Add files file= 'C:\temp2a.sav'
 /file *
 /by surname initial sex dobirth tooth quadrant fillrec.

* Remove additional fee records.

If (missing(trtcode)) trtcode=0.
select if not (range(trtcode,1431,1431,1422,1425,1541,1541,
   1732,1739,1811,1816,2746,2746,5825,5825)).

save outfile='c:\temp3.sav'.

get file='c:\temp3.sav'.

* Reduce to file of teeth not filled in March 2000 schedule but
* where other teeth in the same position were treated then, but treated in following year,
* supplemented by time to failure and
* replacement treatment.

sort cases by surname initial sex dobirth tooth quadrant fillrec doacc trtcode.

aggregate outfile=*
 /presorted
 /break surname initial sex dobirth tooth quadrant fillrec doacc
 / fpcnum persnum partnum suffix postcode docomp exempbox
   trtcode quadcode toothpat exempcod =
 first(fpcnum persnum partnum suffix postcode docomp exempbox
   trtcode quadcode toothpat exempcod)
 /t1=first(trtcode)
 /t2=sum(trtcode)
 /t3=last(trtcode)
 /treats=n.

EXECUTE.

* Calculate time lag for each replacement treatment - up to three per patient.
Do if (fillrec=0 and lag(fillrec)=-1).
compute time=ctime.days(doacc-lag(docomp)).
compute reptrt1=t1.
if (treats=2) reptrt2=t3.
if (treats=3) reptrt2=t2-t1-t3.
if (treats=3) reptrt3=t3.
Compute reptreat=reptrt1.
If (treats>1) reptreat=9999.
End if.

Do if (fillrec=0 and lag(fillrec,2)=-1).
compute time=ctime.days(doacc-lag(docomp,2)).
compute reptrt1=t1.
if (treats=2) reptrt2=t3.
if (treats=3) reptrt2=t2-t1-t3.
if (treats=3) reptrt3=t3.
Compute reptreat=reptrt1.
If (treats>1) reptreat=9999.
End if.

Do if (fillrec=0 and not(lag(fillrec))=-1 and lag(fillrec,3)=-1).
compute time=ctime.days(doacc-lag(docomp,3)).
compute reptrt1=t1.
if (treats=2) reptrt2=t3.
if (treats=3) reptrt2=t2-t1-t3.
if (treats=3) reptrt3=t3.
Compute reptreat=reptrt1.
If (treats>1) reptreat=9999.
End if.

formats trtcode reptreat (f4.0).
formats trtcode reptreat (f4.0) tooth (f1.0).

EXECUTE.

* Tidy up the remaining anomalies.
* Assume negative time cases reflect claims out of sequence.
* Censor all cases at 365 days.
Compute failed=0.
Do if (range(time,1,365)).
compute failed=1.
else.
compute time=365.
end if.
If (missing(time))time=365.

EXECUTE.

* Reduce to one record, for the first retreatment,
* for each quadrant for which the tooth was not treated in March 2000,
* including a dummy quadrant record if the tooth was not retreated.

sort cases by surname initial sex dobirth tooth quadrant time.

Aggregate outfile=* /presorted
 /break surname initial sex dobirth tooth quadrant
 / doacc fpcnum persnum partnum suffix postcode docomp
   exempbox trtcode exempcod time reptreat reptrt1 reptrt2 reptrt3 failed
 =first(doacc fpcnum persnum partnum suffix postcode docomp
   exempbox trtcode exempcod time reptreat reptrt1 reptrt2 reptrt3 failed).

save outfile='D:\pskl\RMatch1975y2001control.sav'.

Get file='D:\pskl\RMatch1975y2001control.sav'.


* Produce some tables.

* a) Proportion of early failures, by tooth position, quadrant, type of filling, and patient sex.

Compute failedPC=100*failed.

autorecode VARIABLES= quadrant sex
 /INTO quadn sexn.

Value Labels
 Tooth TRTCODE
1 '1'
2 '2'
3 '3'
4 '4'
5 '5'
6 '6'
7 '7'
8 '8'
1401 '1401'
1402 '1402'
1403 '1403'
1404 '1404'
1421 '1421'
1426 '1426'.

Variable Labels Tooth 'Tooth Position'
 /quadrant quadn 'Quadrant'
 /trtcode 'Treatment Code'.

Temporary.
select if (any(quadrant,'LL','AA')).

* Basic Tables.
TABLES
  /FORMAT BLANK MISSING('.')
  /OBSERVATION failedpc
  /FTOTAL $t 'Summary'
  /TABLES (sexn > failedpc +  $t )
  BY (tooth +  $t )
  BY ((STATISTICS))
  /STATISTICS
  mean( (PCT7.2 ))
  count( ( COMMA7.0 )) /TITLE  'Retreatment rates by tooth'+
 ' position and patient sex - LL Quadrant'.

Tables

**Notes**
Output Created		17-FEB-2003 18:17:25
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16692
Syntax		TABLES /FORMAT BLANK MISSING('.') /OBSERVATION failedpc /FTOTAL $t 'Summary' /TABLES (sexn > failedpc + $t ) BY (tooth + $t ) BY ((STATISTICS)) /STATISTICS mean( (PCT7.2 )) count( ( COMMA7.0 )) /TITLE 'Retreatment rates by tooth'+ ' position and patient sex - LL Quadrant'.
Resources	Elapsed Time	0:00:00.81

**Retreatment rates by tooth position and patient sex - LL Quadrant**
Mean
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	.66%	.60%	.60%	.99%	1.71%	3.57%	2.68%	1.49%	2.29%
M	.36%	.19%	.85%	.66%	1.85%	2.10%	2.42%	1.49%	1.67%
Summary	.52%	.40%	.73%	.84%	1.78%	2.93%	2.56%	1.49%	2.00%

**Retreatment rates by tooth position and patient sex - LL Quadrant**
Count
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	603	496	335	704	1,226	2,828	2,127	603	8,922
M	554	514	354	602	1,026	2,191	1,858	671	7,770
Summary	1,157	1,010	689	1,306	2,252	5,019	3,985	1,274	16,692

Temporary.
select if (any(quadrant,'LR','AA')).

* Basic Tables.
TABLES
  /FORMAT BLANK MISSING('.')
  /OBSERVATION failedpc
  /FTOTAL $t 'Summary'
  /TABLES (sexn > failedpc +  $t )
  BY (tooth +  $t )
  BY ((STATISTICS))
  /STATISTICS
  mean( (PCT7.2 ))
  count( ( COMMA7.0 )) /TITLE  'Retreatment rates by tooth'+
 ' position and patient sex - LR Quadrant'.

Tables

**Notes**
Output Created		17-FEB-2003 18:17:25
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16704
Syntax		TABLES /FORMAT BLANK MISSING('.') /OBSERVATION failedpc /FTOTAL $t 'Summary' /TABLES (sexn > failedpc + $t ) BY (tooth + $t ) BY ((STATISTICS)) /STATISTICS mean( (PCT7.2 )) count( ( COMMA7.0 )) /TITLE 'Retreatment rates by tooth'+ ' position and patient sex - LR Quadrant'.
Resources	Elapsed Time	0:00:00.60

**Retreatment rates by tooth position and patient sex - LR Quadrant**
Mean
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	.50%	.80%	1.48%	.85%	3.14%	3.25%	2.54%	1.66%	2.38%
M	.72%	1.35%	1.68%	.50%	1.38%	2.01%	2.58%	1.79%	1.78%
Summary	.60%	1.08%	1.58%	.69%	2.35%	2.71%	2.56%	1.73%	2.10%

**Retreatment rates by tooth position and patient sex - LR Quadrant**
Count
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	602	498	337	705	1,241	2,828	2,129	602	8,942
M	558	518	358	597	1,015	2,187	1,857	672	7,762
Summary	1,160	1,016	695	1,302	2,256	5,015	3,986	1,274	16,704

.

Temporary.
select if (any(quadrant,'UL','AA')).

* Basic Tables.
TABLES
  /FORMAT BLANK MISSING('.')
  /OBSERVATION failedpc
  /FTOTAL $t 'Summary'
  /TABLES (sexn > failedpc +  $t )
  BY (tooth +  $t )
  BY ((STATISTICS))
  /STATISTICS
  mean( (PCT7.2 ))
  count( ( COMMA7.0 )) /TITLE  'Retreatment rates by tooth'+
 ' position and patient sex - UL Quadrant'.

Tables

**Notes**
Output Created		17-FEB-2003 18:17:26
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16751
Syntax		TABLES /FORMAT BLANK MISSING('.') /OBSERVATION failedpc /FTOTAL $t 'Summary' /TABLES (sexn > failedpc + $t ) BY (tooth + $t ) BY ((STATISTICS)) /STATISTICS mean( (PCT7.2 )) count( ( COMMA7.0 )) /TITLE 'Retreatment rates by tooth'+ ' position and patient sex - UL Quadrant'.
Resources	Elapsed Time	0:00:00.59

**Retreatment rates by tooth position and patient sex - UL Quadrant**
Mean
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	4.14%	2.37%	2.36%	2.24%	1.96%	2.47%	2.27%	2.31%	2.43%
M	1.95%	3.93%	3.02%	1.80%	1.37%	2.86%	1.90%	2.08%	2.30%
Summary	3.11%	3.17%	2.70%	2.04%	1.69%	2.64%	2.09%	2.19%	2.37%

**Retreatment rates by tooth position and patient sex - UL Quadrant**
Count
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	628	507	339	714	1,225	2,796	2,118	607	8,934
M	563	535	364	612	1,021	2,204	1,845	673	7,817
Summary	1,191	1,042	703	1,326	2,246	5,000	3,963	1,280	16,751

Temporary.
select if (any(quadrant,'UR','AA')).

* Basic Tables.
TABLES
  /FORMAT BLANK MISSING('.')
  /OBSERVATION failedpc
  /FTOTAL $t 'Summary'
  /TABLES (sexn > failedpc +  $t )
  BY (tooth +  $t )
  BY ((STATISTICS))
  /STATISTICS
  mean( (PCT7.2 ))
  count( ( COMMA7.0 )) /TITLE  'Retreatment rates by tooth'+
 ' position and patient sex - UR Quadrant'.

Tables

**Notes**
Output Created		17-FEB-2003 18:17:27
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16699
Syntax		TABLES /FORMAT BLANK MISSING('.') /OBSERVATION failedpc /FTOTAL $t 'Summary' /TABLES (sexn > failedpc + $t ) BY (tooth + $t ) BY ((STATISTICS)) /STATISTICS mean( (PCT7.2 )) count( ( COMMA7.0 )) /TITLE 'Retreatment rates by tooth'+ ' position and patient sex - UR Quadrant'.
Resources	Elapsed Time	0:00:01.09

**Retreatment rates by tooth position and patient sex - UR Quadrant**
Mean
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	2.27%	2.16%	2.08%	2.65%	1.80%	2.71%	2.46%	2.62%	2.43%
M	1.77%	1.92%	2.49%	1.32%	1.77%	1.79%	2.22%	1.05%	1.83%
Summary	2.03%	2.03%	2.29%	2.05%	1.78%	2.31%	2.35%	1.80%	2.15%

**Retreatment rates by tooth position and patient sex - UR Quadrant**
Count
	Tooth Position								Summary
	1	2	3	4	5	6	7	8	Summary
F	616	510	336	716	1,224	2,806	2,116	610	8,934
M	565	522	362	604	1,019	2,178	1,848	667	7,765
Summary	1,181	1,032	698	1,320	2,243	4,984	3,964	1,277	16,699

* b) Replacement treatment , by quadrant.

string retcomb (a14).
compute retcomb=concat(string( reptrt1,f4),'/',string( reptrt2,f4),'/',string( reptrt3,f4)).

Variable Labels retcomb 'Retreatment combination'.

execute.

temporary.
select if ( failed=1).

CROSSTABS
  /TABLES=retcomb  BY quadrant  BY tooth
  /FORMAT= AVALUE TABLES
  /CELLS= COUNT COLUMN .

Crosstabs

**Notes**
Output Created		17-FEB-2003 18:17:29
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	1441
Missing Value Handling	Definition of Missing	User-defined missing values are treated as missing.
Missing Value Handling	Cases Used	Statistics for each table are based on all the cases with valid data in the specified range(s) for all variables in each table.
Syntax		CROSSTABS /TABLES=retcomb BY quadrant BY tooth /FORMAT= AVALUE TABLES /CELLS= COUNT COLUMN .
Resources	Dimensions Requested	3
	Cells Available	80659
	Elapsed Time	0:00:00.81

**Case Processing Summary**
	Cases
	Valid		Missing		Total
	N	Percent	N	Percent	N	Percent
Retreatment combination * quadn 'Quadrant' * Tooth Position	1441	100.0%	0	.0%	1441	100.0%

**Retreatment combination * quadn 'Quadrant' * Tooth Position Crosstabulation**
				quadn 'Quadrant'				Total
Tooth Position				LL	LR	UL	UR	Total
1	Retreatment combination	1401/1501/ .	Count		1			1
		1401/1501/ .	% within quadn 'Quadrant'		14.3%			1.4%
		1421/ ./ .	Count	3	3	23	12	41
		1421/ ./ .	% within quadn 'Quadrant'	50.0%	42.9%	62.2%	50.0%	55.4%
		1421/1501/ .	Count		1	1		2
		1421/1501/ .	% within quadn 'Quadrant'		14.3%	2.7%		2.7%
		1421/1831/ .	Count			1		1
		1421/1831/ .	% within quadn 'Quadrant'			2.7%		1.4%
		1421/5824/ .	Count			1		1
		1421/5824/ .	% within quadn 'Quadrant'			2.7%		1.4%
		1426/ ./ .	Count	1			1	2
		1426/ ./ .	% within quadn 'Quadrant'	16.7%			4.2%	2.7%
		1501/ ./ .	Count			1		1
		1501/ ./ .	% within quadn 'Quadrant'			2.7%		1.4%
		1521/ ./ .	Count			1		1
		1521/ ./ .	% within quadn 'Quadrant'			2.7%		1.4%
		1601/ ./ .	Count				5	5
		1601/ ./ .	% within quadn 'Quadrant'				20.8%	6.8%
		1601/3611/ .	Count			1		1
		1601/3611/ .	% within quadn 'Quadrant'			2.7%		1.4%
		1721/ ./ .	Count			2	1	3
		1721/ ./ .	% within quadn 'Quadrant'			5.4%	4.2%	4.1%
		1782/ ./ .	Count				1	1
		1782/ ./ .	% within quadn 'Quadrant'				4.2%	1.4%
		1825/ ./ .	Count			1		1
		1825/ ./ .	% within quadn 'Quadrant'			2.7%		1.4%
		1826/ ./ .	Count				1	1
		1826/ ./ .	% within quadn 'Quadrant'				4.2%	1.4%
		2101/ ./ .	Count			1	1	2
		2101/ ./ .	% within quadn 'Quadrant'			2.7%	4.2%	2.7%
		2733/ ./ .	Count	1	1	1	1	4
		2733/ ./ .	% within quadn 'Quadrant'	16.7%	14.3%	2.7%	4.2%	5.4%
		3611/ ./ .	Count	1	1	3	1	6
		3611/ ./ .	% within quadn 'Quadrant'	16.7%	14.3%	8.1%	4.2%	8.1%
	Total		Count	6	7	37	24	74
	Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%
2	Retreatment combination	1421/ ./ .	Count	1	4	18	7	30
		1421/ ./ .	% within quadn 'Quadrant'	25.0%	36.4%	54.5%	33.3%	43.5%
		1421/1501/ .	Count			1		1
		1421/1501/ .	% within quadn 'Quadrant'			3.0%		1.4%
		1426/ ./ .	Count	2	3	3	1	9
		1426/ ./ .	% within quadn 'Quadrant'	50.0%	27.3%	9.1%	4.8%	13.0%
		1426/1501/ .	Count				1	1
		1426/1501/ .	% within quadn 'Quadrant'				4.8%	1.4%
		1501/1721/ .	Count			1		1
		1501/1721/ .	% within quadn 'Quadrant'			3.0%		1.4%
		1721/ ./ .	Count		1	1	4	6
		1721/ ./ .	% within quadn 'Quadrant'		9.1%	3.0%	19.0%	8.7%
		1782/ ./ .	Count		1	2	5	8
		1782/ ./ .	% within quadn 'Quadrant'		9.1%	6.1%	23.8%	11.6%
		1832/ ./ .	Count			1		1
		1832/ ./ .	% within quadn 'Quadrant'			3.0%		1.4%
		2101/ ./ .	Count			1		1
		2101/ ./ .	% within quadn 'Quadrant'			3.0%		1.4%
		2733/ ./ .	Count	1	2	3	2	8
		2733/ ./ .	% within quadn 'Quadrant'	25.0%	18.2%	9.1%	9.5%	11.6%
		3611/ ./ .	Count			2	1	3
		3611/ ./ .	% within quadn 'Quadrant'			6.1%	4.8%	4.3%
	Total		Count	4	11	33	21	69
	Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%
3	Retreatment combination	1421/ ./ .	Count	4	4	9	9	26
		1421/ ./ .	% within quadn 'Quadrant'	80.0%	36.4%	47.4%	56.3%	51.0%
		1421/1501/ .	Count				2	2
		1421/1501/ .	% within quadn 'Quadrant'				12.5%	3.9%
		1426/ ./ .	Count	1	3	2	2	8
		1426/ ./ .	% within quadn 'Quadrant'	20.0%	27.3%	10.5%	12.5%	15.7%
		1721/ ./ .	Count			1	1	2
		1721/ ./ .	% within quadn 'Quadrant'			5.3%	6.3%	3.9%
		1782/ ./ .	Count			1		1
		1782/ ./ .	% within quadn 'Quadrant'			5.3%		2.0%
		1807/ ./ .	Count			1		1
		1807/ ./ .	% within quadn 'Quadrant'			5.3%		2.0%
		2101/ ./ .	Count		1	2		3
		2101/ ./ .	% within quadn 'Quadrant'		9.1%	10.5%		5.9%
		2733/ ./ .	Count		1	3	2	6
		2733/ ./ .	% within quadn 'Quadrant'		9.1%	15.8%	12.5%	11.8%
		2863/ ./ .	Count		1			1
		2863/ ./ .	% within quadn 'Quadrant'		9.1%			2.0%
		6331/ ./ .	Count		1			1
		6331/ ./ .	% within quadn 'Quadrant'		9.1%			2.0%
	Total		Count	5	11	19	16	51
	Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%
4	Retreatment combination	1401/ ./ .	Count				1	1
		1401/ ./ .	% within quadn 'Quadrant'				3.7%	1.4%
		1403/ ./ .	Count	2	3	12	9	26
		1403/ ./ .	% within quadn 'Quadrant'	18.2%	33.3%	44.4%	33.3%	35.1%
		1403/1426/ .	Count			1		1
		1403/1426/ .	% within quadn 'Quadrant'			3.7%		1.4%
		1403/1501/ .	Count	1			1	2
		1403/1501/ .	% within quadn 'Quadrant'	9.1%			3.7%	2.7%
		1404/ ./ .	Count			2	4	6
		1404/ ./ .	% within quadn 'Quadrant'			7.4%	14.8%	8.1%
		1411/ ./ .	Count			1		1
		1411/ ./ .	% within quadn 'Quadrant'			3.7%		1.4%
		1421/ ./ .	Count	5	1	4	6	16
		1421/ ./ .	% within quadn 'Quadrant'	45.5%	11.1%	14.8%	22.2%	21.6%
		1421/1426/ .	Count		1			1
		1421/1426/ .	% within quadn 'Quadrant'		11.1%			1.4%
		1426/ ./ .	Count	1	1	1		3
		1426/ ./ .	% within quadn 'Quadrant'	9.1%	11.1%	3.7%		4.1%
		1441/ ./ .	Count	1		1	1	3
		1441/ ./ .	% within quadn 'Quadrant'	9.1%		3.7%	3.7%	4.1%
		1441/1421/ .	Count				1	1
		1441/1421/ .	% within quadn 'Quadrant'				3.7%	1.4%
		1442/ ./ .	Count				1	1
		1442/ ./ .	% within quadn 'Quadrant'				3.7%	1.4%
		1443/ ./ .	Count		2			2
		1443/ ./ .	% within quadn 'Quadrant'		22.2%			2.7%
		1444/ ./ .	Count	1		1		2
		1444/ ./ .	% within quadn 'Quadrant'	9.1%		3.7%		2.7%
		1704/ ./ .	Count				1	1
		1704/ ./ .	% within quadn 'Quadrant'				3.7%	1.4%
		1721/ ./ .	Count			1		1
		1721/ ./ .	% within quadn 'Quadrant'			3.7%		1.4%
		1782/ ./ .	Count			1		1
		1782/ ./ .	% within quadn 'Quadrant'			3.7%		1.4%
		1826/ ./ .	Count			1		1
		1826/ ./ .	% within quadn 'Quadrant'			3.7%		1.4%
		2733/ ./ .	Count		1	1	1	3
		2733/ ./ .	% within quadn 'Quadrant'		11.1%	3.7%	3.7%	4.1%
		2745/ ./ .	Count				1	1
		2745/ ./ .	% within quadn 'Quadrant'				3.7%	1.4%
	Total		Count	11	9	27	27	74
	Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%
5	Retreatment combination	1401/ ./ .	Count	2		2	1	5
		1401/ ./ .	% within quadn 'Quadrant'	5.0%		5.3%	2.5%	2.9%
		1401/1421/ .	Count	1				1
		1401/1421/ .	% within quadn 'Quadrant'	2.5%				.6%
		1401/1426/ .	Count		1			1
		1401/1426/ .	% within quadn 'Quadrant'		1.9%			.6%
		1401/1501/ .	Count			1		1
		1401/1501/ .	% within quadn 'Quadrant'			2.6%		.6%
		1402/ ./ .	Count			1		1
		1402/ ./ .	% within quadn 'Quadrant'			2.6%		.6%
		1403/ ./ .	Count	12	18	12	13	55
		1403/ ./ .	% within quadn 'Quadrant'	30.0%	34.0%	31.6%	32.5%	32.2%
		1403/1426/ .	Count		1			1
		1403/1426/ .	% within quadn 'Quadrant'		1.9%			.6%
		1403/1501/ .	Count	2	2	2	4	10
		1403/1501/ .	% within quadn 'Quadrant'	5.0%	3.8%	5.3%	10.0%	5.8%
		1404/ ./ .	Count	5	2	6	4	17
		1404/ ./ .	% within quadn 'Quadrant'	12.5%	3.8%	15.8%	10.0%	9.9%
		1404/1426/ .	Count			1		1
		1404/1426/ .	% within quadn 'Quadrant'			2.6%		.6%
		1404/1501/ .	Count	2			1	3
		1404/1501/ .	% within quadn 'Quadrant'	5.0%			2.5%	1.8%
		1421/ ./ .	Count	3	2		4	9
		1421/ ./ .	% within quadn 'Quadrant'	7.5%	3.8%		10.0%	5.3%
		1421/1501/ .	Count				1	1
		1421/1501/ .	% within quadn 'Quadrant'				2.5%	.6%
		1426/ ./ .	Count	3	5	1	2	11
		1426/ ./ .	% within quadn 'Quadrant'	7.5%	9.4%	2.6%	5.0%	6.4%
		1441/ ./ .	Count			1		1
		1441/ ./ .	% within quadn 'Quadrant'			2.6%		.6%
		1443/ ./ .	Count		1		1	2
		1443/ ./ .	% within quadn 'Quadrant'		1.9%		2.5%	1.2%
		1444/ ./ .	Count		1	1		2
		1444/ ./ .	% within quadn 'Quadrant'		1.9%	2.6%		1.2%
		1501/ ./ .	Count		1			1
		1501/ ./ .	% within quadn 'Quadrant'		1.9%			.6%
		1701/ ./ .	Count			1	1	2
		1701/ ./ .	% within quadn 'Quadrant'			2.6%	2.5%	1.2%
		1704/ ./ .	Count				1	1
		1704/ ./ .	% within quadn 'Quadrant'				2.5%	.6%
		1721/ ./ .	Count	3	4	2	2	11
		1721/ ./ .	% within quadn 'Quadrant'	7.5%	7.5%	5.3%	5.0%	6.4%
		1782/ ./ .	Count	1	2			3
		1782/ ./ .	% within quadn 'Quadrant'	2.5%	3.8%			1.8%
		1782/1744/ .	Count	1				1
		1782/1744/ .	% within quadn 'Quadrant'	2.5%				.6%
		1807/ ./ .	Count		2			2
		1807/ ./ .	% within quadn 'Quadrant'		3.8%			1.2%
		2101/ ./ .	Count	2	2	2	1	7
		2101/ ./ .	% within quadn 'Quadrant'	5.0%	3.8%	5.3%	2.5%	4.1%
		2201/ ./ .	Count		1		2	3
		2201/ ./ .	% within quadn 'Quadrant'		1.9%		5.0%	1.8%
		2203/ ./ .	Count	1		1	1	3
		2203/ ./ .	% within quadn 'Quadrant'	2.5%		2.6%	2.5%	1.8%
		2733/ ./ .	Count		3	3	1	7
		2733/ ./ .	% within quadn 'Quadrant'		5.7%	7.9%	2.5%	4.1%
		2744/ ./ .	Count	1	1			2
		2744/ ./ .	% within quadn 'Quadrant'	2.5%	1.9%			1.2%
		3611/ ./ .	Count		1			1
		3611/ ./ .	% within quadn 'Quadrant'		1.9%			.6%
		3611/1704/ .	Count		1			1
		3611/1704/ .	% within quadn 'Quadrant'		1.9%			.6%
		5813/ ./ .	Count			1		1
		5813/ ./ .	% within quadn 'Quadrant'			2.6%		.6%
6301/ ./ .	Count		1			1
6301/ ./ .	% within quadn 'Quadrant'		1.9%			.6%
6341/ ./ .	Count	1	1			2
6341/ ./ .	% within quadn 'Quadrant'	2.5%	1.9%			1.2%
Total		Count	40	53	38	40	171
Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%
6	Retreatment combination	1131/1711/ .	Count		1			1
		1131/1711/ .	% within quadn 'Quadrant'		.7%			.2%
		1401/ ./ .	Count	21	20	22	21	84
		1401/ ./ .	% within quadn 'Quadrant'	14.3%	14.7%	16.7%	18.3%	15.8%
		1401/1403/ .	Count			1		1
		1401/1403/ .	% within quadn 'Quadrant'			.8%		.2%
		1401/1421/ .	Count	1	1		1	3
		1401/1421/ .	% within quadn 'Quadrant'	.7%	.7%		.9%	.6%
		1401/1426/ .	Count		1			1
		1401/1426/ .	% within quadn 'Quadrant'		.7%			.2%
		1401/1501/ .	Count	1				1
		1401/1501/ .	% within quadn 'Quadrant'	.7%				.2%
		1402/ ./ .	Count	10	3	4	5	22
		1402/ ./ .	% within quadn 'Quadrant'	6.8%	2.2%	3.0%	4.3%	4.2%
		1402/1501/ .	Count	1	1			2
		1402/1501/ .	% within quadn 'Quadrant'	.7%	.7%			.4%
		1403/ ./ .	Count	39	39	52	40	170
		1403/ ./ .	% within quadn 'Quadrant'	26.5%	28.7%	39.4%	34.8%	32.1%
		1403/1421/ .	Count	1	1			2
		1403/1421/ .	% within quadn 'Quadrant'	.7%	.7%			.4%
		1403/1426/ .	Count	1				1
		1403/1426/ .	% within quadn 'Quadrant'	.7%				.2%
		1403/1501/ .	Count	3	5	3	4	15
		1403/1501/ .	% within quadn 'Quadrant'	2.0%	3.7%	2.3%	3.5%	2.8%
		1403/1825/ .	Count				1	1
		1403/1825/ .	% within quadn 'Quadrant'				.9%	.2%
		1404/ ./ .	Count	14	12	8	12	46
		1404/ ./ .	% within quadn 'Quadrant'	9.5%	8.8%	6.1%	10.4%	8.7%
		1404/1501/ .	Count	1	2	1	2	6
		1404/1501/ .	% within quadn 'Quadrant'	.7%	1.5%	.8%	1.7%	1.1%
		1411/ ./ .	Count	1		1		2
		1411/ ./ .	% within quadn 'Quadrant'	.7%		.8%		.4%
		1421/ ./ .	Count	6	7	5	5	23
		1421/ ./ .	% within quadn 'Quadrant'	4.1%	5.1%	3.8%	4.3%	4.3%
		1421/1442/ .	Count			1		1
		1421/1442/ .	% within quadn 'Quadrant'			.8%		.2%
		1426/ ./ .	Count	6	8	8	3	25
		1426/ ./ .	% within quadn 'Quadrant'	4.1%	5.9%	6.1%	2.6%	4.7%
		1426/1444/ .	Count	1				1
		1426/1444/ .	% within quadn 'Quadrant'	.7%				.2%
		1426/1501/ .	Count		1			1
		1426/1501/ .	% within quadn 'Quadrant'		.7%			.2%
		1441/ ./ .	Count	3	1	2	2	8
		1441/ ./ .	% within quadn 'Quadrant'	2.0%	.7%	1.5%	1.7%	1.5%
		1442/ ./ .	Count				1	1
		1442/ ./ .	% within quadn 'Quadrant'				.9%	.2%
		1443/ ./ .	Count	1		1	2	4
		1443/ ./ .	% within quadn 'Quadrant'	.7%		.8%	1.7%	.8%
		1444/ ./ .	Count	1	2	2	3	8
		1444/ ./ .	% within quadn 'Quadrant'	.7%	1.5%	1.5%	2.6%	1.5%
		1501/ ./ .	Count		1	3		4
		1501/ ./ .	% within quadn 'Quadrant'		.7%	2.3%		.8%
		1501/1711/ .	Count	2				2
		1501/1711/ .	% within quadn 'Quadrant'	1.4%				.4%
		1702/ ./ .	Count		1			1
		1702/ ./ .	% within quadn 'Quadrant'		.7%			.2%
		1704/ ./ .	Count	4	2	2	1	9
		1704/ ./ .	% within quadn 'Quadrant'	2.7%	1.5%	1.5%	.9%	1.7%
		1711/ ./ .	Count	2	3	1	2	8
		1711/ ./ .	% within quadn 'Quadrant'	1.4%	2.2%	.8%	1.7%	1.5%
		1781/ ./ .	Count	1	1			2
		1781/ ./ .	% within quadn 'Quadrant'	.7%	.7%			.4%
		1782/ ./ .	Count		2		1	3
		1782/ ./ .	% within quadn 'Quadrant'		1.5%		.9%	.6%
1825/ ./ .	Count		2			2
1825/ ./ .	% within quadn 'Quadrant'		1.5%			.4%
1826/ ./ .	Count		1			1
1826/ ./ .	% within quadn 'Quadrant'		.7%			.2%
2101/ ./ .	Count	11	3	6	2	22
2101/ ./ .	% within quadn 'Quadrant'	7.5%	2.2%	4.5%	1.7%	4.2%
2201/ ./ .	Count	2	1			3
2201/ ./ .	% within quadn 'Quadrant'	1.4%	.7%			.6%
2203/ ./ .	Count	4	6	1	1	12
2203/ ./ .	% within quadn 'Quadrant'	2.7%	4.4%	.8%	.9%	2.3%
2733/ ./ .	Count	1	2	4	1	8
2733/ ./ .	% within quadn 'Quadrant'	.7%	1.5%	3.0%	.9%	1.5%
3611/ ./ .	Count	4	2	2	4	12
3611/ ./ .	% within quadn 'Quadrant'	2.7%	1.5%	1.5%	3.5%	2.3%
5001/ ./ .	Count	1	1	1		3
5001/ ./ .	% within quadn 'Quadrant'	.7%	.7%	.8%		.6%
5201/ ./ .	Count	1				1
5201/ ./ .	% within quadn 'Quadrant'	.7%				.2%
6264/ ./ .	Count		1			1
6264/ ./ .	% within quadn 'Quadrant'		.7%			.2%
6301/ ./ .	Count	2	1	1		4
6301/ ./ .	% within quadn 'Quadrant'	1.4%	.7%	.8%		.8%
6332/ ./ .	Count				1	1
6332/ ./ .	% within quadn 'Quadrant'				.9%	.2%
6341/ ./ .	Count		1			1
6341/ ./ .	% within quadn 'Quadrant'		.7%			.2%
Total		Count	147	136	132	115	530
Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%
7	Retreatment combination	1401/ ./ .	Count	31	31	30	21	113
		1401/ ./ .	% within quadn 'Quadrant'	30.4%	30.4%	36.1%	22.6%	29.7%
		1401/1421/ .	Count		1	1		2
		1401/1421/ .	% within quadn 'Quadrant'		1.0%	1.2%		.5%
		1401/1426/ .	Count		1			1
		1401/1426/ .	% within quadn 'Quadrant'		1.0%			.3%
		1401/1501/ .	Count			1		1
		1401/1501/ .	% within quadn 'Quadrant'			1.2%		.3%
		1401/1501/6331	Count		1			1
		1401/1501/6331	% within quadn 'Quadrant'		1.0%			.3%
		1402/ ./ .	Count	4	4	2	7	17
		1402/ ./ .	% within quadn 'Quadrant'	3.9%	3.9%	2.4%	7.5%	4.5%
		1403/ ./ .	Count	23	20	22	28	93
		1403/ ./ .	% within quadn 'Quadrant'	22.5%	19.6%	26.5%	30.1%	24.5%
		1403/1421/ .	Count				1	1
		1403/1421/ .	% within quadn 'Quadrant'				1.1%	.3%
		1403/1426/ .	Count			1		1
		1403/1426/ .	% within quadn 'Quadrant'			1.2%		.3%
		1403/1501/ .	Count	1		1		2
		1403/1501/ .	% within quadn 'Quadrant'	1.0%		1.2%		.5%
		1404/ ./ .	Count	6	2	5	4	17
		1404/ ./ .	% within quadn 'Quadrant'	5.9%	2.0%	6.0%	4.3%	4.5%
		1404/1501/ .	Count		1		1	2
		1404/1501/ .	% within quadn 'Quadrant'		1.0%		1.1%	.5%
		1421/ ./ .	Count	6	4	2	3	15
		1421/ ./ .	% within quadn 'Quadrant'	5.9%	3.9%	2.4%	3.2%	3.9%
		1426/ ./ .	Count	6	14	4	4	28
		1426/ ./ .	% within quadn 'Quadrant'	5.9%	13.7%	4.8%	4.3%	7.4%
		1426/1501/ .	Count				1	1
		1426/1501/ .	% within quadn 'Quadrant'				1.1%	.3%
		1441/ ./ .	Count	1	2	2	5	10
		1441/ ./ .	% within quadn 'Quadrant'	1.0%	2.0%	2.4%	5.4%	2.6%
		1442/ ./ .	Count	1	2		1	4
		1442/ ./ .	% within quadn 'Quadrant'	1.0%	2.0%		1.1%	1.1%
		1443/ ./ .	Count	1	1		1	3
		1443/ ./ .	% within quadn 'Quadrant'	1.0%	1.0%		1.1%	.8%
		1444/ ./ .	Count	4	2	3	5	14
		1444/ ./ .	% within quadn 'Quadrant'	3.9%	2.0%	3.6%	5.4%	3.7%
		1501/ ./ .	Count	1		1		2
		1501/ ./ .	% within quadn 'Quadrant'	1.0%		1.2%		.5%
		1702/ ./ .	Count	1				1
		1702/ ./ .	% within quadn 'Quadrant'	1.0%				.3%
		1704/ ./ .	Count	2				2
		1704/ ./ .	% within quadn 'Quadrant'	2.0%				.5%
		1711/ ./ .	Count	2	2	1	1	6
		1711/ ./ .	% within quadn 'Quadrant'	2.0%	2.0%	1.2%	1.1%	1.6%
		1807/ ./ .	Count		1			1
		1807/ ./ .	% within quadn 'Quadrant'		1.0%			.3%
		2101/ ./ .	Count	3	5		2	10
		2101/ ./ .	% within quadn 'Quadrant'	2.9%	4.9%		2.2%	2.6%
		2201/ ./ .	Count			2	2	4
		2201/ ./ .	% within quadn 'Quadrant'			2.4%	2.2%	1.1%
		2203/ ./ .	Count	2	4	1	3	10
		2203/ ./ .	% within quadn 'Quadrant'	2.0%	3.9%	1.2%	3.2%	2.6%
		2733/ ./ .	Count	1		2		3
		2733/ ./ .	% within quadn 'Quadrant'	1.0%		2.4%		.8%
		3611/ ./ .	Count	3	1	1	1	6
		3611/ ./ .	% within quadn 'Quadrant'	2.9%	1.0%	1.2%	1.1%	1.6%
		3611/1501/ .	Count				1	1
		3611/1501/ .	% within quadn 'Quadrant'				1.1%	.3%
		6301/ ./ .	Count	2	1	1	1	5
		6301/ ./ .	% within quadn 'Quadrant'	2.0%	1.0%	1.2%	1.1%	1.3%
		6332/ ./ .	Count	1	1			2
		6332/ ./ .	% within quadn 'Quadrant'	1.0%	1.0%			.5%
6341/1501/ .	Count		1			1
6341/1501/ .	% within quadn 'Quadrant'		1.0%			.3%
Total		Count	102	102	83	93	380
Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%
8	Retreatment combination	1401/ ./ .	Count	7	9	10	8	34
		1401/ ./ .	% within quadn 'Quadrant'	36.8%	40.9%	35.7%	34.8%	37.0%
		1401/1426/ .	Count		1			1
		1401/1426/ .	% within quadn 'Quadrant'		4.5%			1.1%
		1402/ ./ .	Count	1	2	1	1	5
		1402/ ./ .	% within quadn 'Quadrant'	5.3%	9.1%	3.6%	4.3%	5.4%
		1403/ ./ .	Count	4	2	4	3	13
		1403/ ./ .	% within quadn 'Quadrant'	21.1%	9.1%	14.3%	13.0%	14.1%
		1404/ ./ .	Count	1	2			3
		1404/ ./ .	% within quadn 'Quadrant'	5.3%	9.1%			3.3%
		1426/ ./ .	Count	1	1	1	2	5
		1426/ ./ .	% within quadn 'Quadrant'	5.3%	4.5%	3.6%	8.7%	5.4%
		1441/ ./ .	Count	1	1			2
		1441/ ./ .	% within quadn 'Quadrant'	5.3%	4.5%			2.2%
		1442/ ./ .	Count		1			1
		1442/ ./ .	% within quadn 'Quadrant'		4.5%			1.1%
		1711/ ./ .	Count	1				1
		1711/ ./ .	% within quadn 'Quadrant'	5.3%				1.1%
		2101/ ./ .	Count	2	1	6	6	15
		2101/ ./ .	% within quadn 'Quadrant'	10.5%	4.5%	21.4%	26.1%	16.3%
		2201/ ./ .	Count			3	2	5
		2201/ ./ .	% within quadn 'Quadrant'			10.7%	8.7%	5.4%
		2203/ ./ .	Count			2	1	3
		2203/ ./ .	% within quadn 'Quadrant'			7.1%	4.3%	3.3%
		2204/ ./ .	Count		2			2
		2204/ ./ .	% within quadn 'Quadrant'		9.1%			2.2%
		2205/ ./ .	Count	1				1
		2205/ ./ .	% within quadn 'Quadrant'	5.3%				1.1%
		6301/ ./ .	Count			1		1
		6301/ ./ .	% within quadn 'Quadrant'			3.6%		1.1%
	Total		Count	19	22	28	23	92
	Total		% within quadn 'Quadrant'	100.0%	100.0%	100.0%	100.0%	100.0%

* c) Kaplan-Meier Survivor Functions and Hazard functions.

temporary.
select if (any(quadrant,'LL','AA')).
KM
  time  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL.

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:30
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16692
Syntax		KM time /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL.
Resources	Elapsed Time	0:00:02.30

temporary.
select if (any(quadrant,'LR','AA')).
KM
  time  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL.

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:32
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16704
Syntax		KM time /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL.
Resources	Elapsed Time	0:00:01.05

temporary.
select if (any(quadrant,'UL','AA')).
KM
  time  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL.

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:33
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16751
Syntax		KM time /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL.
Resources	Elapsed Time	0:00:01.12

temporary.
select if (any(quadrant,'UR','AA')).
KM
  time  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL.

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:35
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16699
Syntax		KM time /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL.
Resources	Elapsed Time	0:00:01.33

* Subdivide by tooth.

Get file='D:\pskl\RMatch1975y2001control.sav'.

select if (any(quadrant,'LL','AA')).
KM
  time  BY tooth  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL
  /save Survival(survtoot).

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:36
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16692
Syntax		KM time BY tooth /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL /save Survival(survtoot).
Resources	Elapsed Time	0:00:01.55

 Survival Analysis for TIME

                            Total        Number         Number          Percent
                                         Events        Censored        Censored

  TOOTH           1          1157             6           1151            99.48
  TOOTH           2          1010             4           1006            99.60
  TOOTH           3           689             5            684            99.27
  TOOTH           4          1306            11           1295            99.16
  TOOTH           5          2252            40           2212            98.22
  TOOTH           6          5019           147           4872            97.07
  TOOTH           7          3985           102           3883            97.44
  TOOTH           8          1274            19           1255            98.51

Overall                     16692           334          16358            98.00

* Create Excel file for plotting.

select if (failed=1).
sort cases by time tooth.

aggregate outfile=*
 /presorted
 /break time tooth
 /survtoot
  =first(survtoot)
 /n=n.

If (tooth=1) tooth1=survtoot.
If (tooth=2) tooth2=survtoot.
If (tooth=3) tooth3=survtoot.
If (tooth=4) tooth4=survtoot.
If (tooth=5) tooth5=survtoot.
If (tooth=6) tooth6=survtoot.
If (tooth=7) tooth7=survtoot.
If (tooth=8) tooth8=survtoot.

aggregate outfile=*
 /presorted
 /break time
 / tooth1 to tooth8
 = first(tooth1 to tooth8).

SAVE TRANSLATE OUTFILE='I:\Research Projects\longevity\early life 1975 control LL.xls'
  /TYPE=XLS /MAP /REPLACE /FIELDNAMES.

Data written to I:\Research Projects\longevity\early life 1975 control LL.xls.
9 variables and 201 cases written to range: SPSS.
Variable: TIME       Type: Number   Width:  8   Dec: 2
Variable: TOOTH1     Type: Number   Width:  8   Dec: 2
Variable: TOOTH2     Type: Number   Width:  8   Dec: 2
Variable: TOOTH3     Type: Number   Width:  8   Dec: 2
Variable: TOOTH4     Type: Number   Width:  8   Dec: 2
Variable: TOOTH5     Type: Number   Width:  8   Dec: 2
Variable: TOOTH6     Type: Number   Width:  8   Dec: 2
Variable: TOOTH7     Type: Number   Width:  8   Dec: 2
Variable: TOOTH8     Type: Number   Width:  8   Dec: 2






Get file='D:\pskl\RMatch1975y2001control.sav'.

select if (any(quadrant,'LR','AA')).
KM
  time  BY tooth  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL
  /save Survival(survtoot).

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:39
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16704
Syntax		KM time BY tooth /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL /save Survival(survtoot).
Resources	Elapsed Time	0:00:01.65

 Survival Analysis for TIME

                            Total        Number         Number          Percent
                                         Events        Censored        Censored

  TOOTH           1          1160             7           1153            99.40
  TOOTH           2          1016            11           1005            98.92
  TOOTH           3           695            11            684            98.42
  TOOTH           4          1302             9           1293            99.31
  TOOTH           5          2256            53           2203            97.65
  TOOTH           6          5015           136           4879            97.29
  TOOTH           7          3986           102           3884            97.44
  TOOTH           8          1274            22           1252            98.27

Overall                     16704           351          16353            97.90

* Create Excel file for plotting.

select if (failed=1).
sort cases by time tooth.

aggregate outfile=*
 /presorted
 /break time tooth
 /survtoot
  =first(survtoot)
 /n=n.

If (tooth=1) tooth1=survtoot.
If (tooth=2) tooth2=survtoot.
If (tooth=3) tooth3=survtoot.
If (tooth=4) tooth4=survtoot.
If (tooth=5) tooth5=survtoot.
If (tooth=6) tooth6=survtoot.
If (tooth=7) tooth7=survtoot.
If (tooth=8) tooth8=survtoot.

aggregate outfile=*
 /presorted
 /break time
 / tooth1 to tooth8
 = first(tooth1 to tooth8).

SAVE TRANSLATE OUTFILE='I:\Research Projects\longevity\early life 1975 control LR.xls'
  /TYPE=XLS /MAP /REPLACE /FIELDNAMES.

Data written to I:\Research Projects\longevity\early life 1975 control LR.xls.
9 variables and 213 cases written to range: SPSS.
Variable: TIME       Type: Number   Width:  8   Dec: 2
Variable: TOOTH1     Type: Number   Width:  8   Dec: 2
Variable: TOOTH2     Type: Number   Width:  8   Dec: 2
Variable: TOOTH3     Type: Number   Width:  8   Dec: 2
Variable: TOOTH4     Type: Number   Width:  8   Dec: 2
Variable: TOOTH5     Type: Number   Width:  8   Dec: 2
Variable: TOOTH6     Type: Number   Width:  8   Dec: 2
Variable: TOOTH7     Type: Number   Width:  8   Dec: 2
Variable: TOOTH8     Type: Number   Width:  8   Dec: 2









Get file='D:\pskl\RMatch1975y2001control.sav'.

select if (any(quadrant,'UL','AA')).
KM
  time  BY tooth  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL
  /save Survival(survtoot).

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:41
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16751
Syntax		KM time BY tooth /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL /save Survival(survtoot).
Resources	Elapsed Time	0:00:01.26

 Survival Analysis for TIME

                            Total        Number         Number          Percent
                                         Events        Censored        Censored

  TOOTH           1          1191            37           1154            96.89
  TOOTH           2          1042            33           1009            96.83
  TOOTH           3           703            19            684            97.30
  TOOTH           4          1326            27           1299            97.96
  TOOTH           5          2246            38           2208            98.31
  TOOTH           6          5000           132           4868            97.36
  TOOTH           7          3963            83           3880            97.91
  TOOTH           8          1280            28           1252            97.81

Overall                     16751           397          16354            97.63

* Create Excel file for plotting.

select if (failed=1).
sort cases by time tooth.

aggregate outfile=*
 /presorted
 /break time tooth
 /survtoot
  =first(survtoot)
 /n=n.

If (tooth=1) tooth1=survtoot.
If (tooth=2) tooth2=survtoot.
If (tooth=3) tooth3=survtoot.
If (tooth=4) tooth4=survtoot.
If (tooth=5) tooth5=survtoot.
If (tooth=6) tooth6=survtoot.
If (tooth=7) tooth7=survtoot.
If (tooth=8) tooth8=survtoot.

aggregate outfile=*
 /presorted
 /break time
 / tooth1 to tooth8
 = first(tooth1 to tooth8).

SAVE TRANSLATE OUTFILE='I:\Research Projects\longevity\early life 1975 control UL.xls'
  /TYPE=XLS /MAP /REPLACE /FIELDNAMES.

Data written to I:\Research Projects\longevity\early life 1975 control UL.xls.
9 variables and 213 cases written to range: SPSS.
Variable: TIME       Type: Number   Width:  8   Dec: 2
Variable: TOOTH1     Type: Number   Width:  8   Dec: 2
Variable: TOOTH2     Type: Number   Width:  8   Dec: 2
Variable: TOOTH3     Type: Number   Width:  8   Dec: 2
Variable: TOOTH4     Type: Number   Width:  8   Dec: 2
Variable: TOOTH5     Type: Number   Width:  8   Dec: 2
Variable: TOOTH6     Type: Number   Width:  8   Dec: 2
Variable: TOOTH7     Type: Number   Width:  8   Dec: 2
Variable: TOOTH8     Type: Number   Width:  8   Dec: 2









Get file='D:\pskl\RMatch1975y2001control.sav'.

select if (any(quadrant,'UR','AA')).
KM
  time  BY tooth  /STATUS=failed(1)
  /PRINT NONE
  /PLOT SURVIVAL
  /save Survival(survtoot).

Kaplan-Meier

**Notes**
Output Created		17-FEB-2003 18:17:43
Comments
Input	Data	D:\pskl\RMatch1975y2001control.sav
	File Label	Aggregated File
	Filter	<none>
	Weight	<none>
	Split File	<none>
	N of Rows in Working Data File	16699
Syntax		KM time BY tooth /STATUS=failed(1) /PRINT NONE /PLOT SURVIVAL /save Survival(survtoot).
Resources	Elapsed Time	0:00:02.74

 Survival Analysis for TIME

                            Total        Number         Number          Percent
                                         Events        Censored        Censored

  TOOTH           1          1181            24           1157            97.97
  TOOTH           2          1032            21           1011            97.97
  TOOTH           3           698            16            682            97.71
  TOOTH           4          1320            27           1293            97.95
  TOOTH           5          2243            40           2203            98.22
  TOOTH           6          4984           115           4869            97.69
  TOOTH           7          3964            93           3871            97.65
  TOOTH           8          1277            23           1254            98.20

Overall                     16699           359          16340            97.85

* Create Excel file for plotting.

select if (failed=1).
sort cases by time tooth.

aggregate outfile=*
 /presorted
 /break time tooth
 /survtoot
  =first(survtoot)
 /n=n.

If (tooth=1) tooth1=survtoot.
If (tooth=2) tooth2=survtoot.
If (tooth=3) tooth3=survtoot.
If (tooth=4) tooth4=survtoot.
If (tooth=5) tooth5=survtoot.
If (tooth=6) tooth6=survtoot.
If (tooth=7) tooth7=survtoot.
If (tooth=8) tooth8=survtoot.

aggregate outfile=*
 /presorted
 /break time
 / tooth1 to tooth8
 = first(tooth1 to tooth8).

SAVE TRANSLATE OUTFILE='I:\Research Projects\longevity\early life 1975 control UR.xls'
  /TYPE=XLS /MAP /REPLACE /FIELDNAMES.

Data written to I:\Research Projects\longevity\early life 1975 control UR.xls.
9 variables and 211 cases written to range: SPSS.
Variable: TIME       Type: Number   Width:  8   Dec: 2
Variable: TOOTH1     Type: Number   Width:  8   Dec: 2
Variable: TOOTH2     Type: Number   Width:  8   Dec: 2
Variable: TOOTH3     Type: Number   Width:  8   Dec: 2
Variable: TOOTH4     Type: Number   Width:  8   Dec: 2
Variable: TOOTH5     Type: Number   Width:  8   Dec: 2
Variable: TOOTH6     Type: Number   Width:  8   Dec: 2
Variable: TOOTH7     Type: Number   Width:  8   Dec: 2
Variable: TOOTH8     Type: Number   Width:  8   Dec: 2